Counterintuitive transitions between crossing energy levels 
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We calculate analytically the probabilities for intuitive and counterintuitive transitions in a three- 
state system, in which two parallel energies are crossed by a third, tilted energy. The state with 
the tilted energy is coupled to the other two states in a chainwise linkage pattern with constant 
couplings of finite duration. The probability for a counterintuitive transition is found to increase 
with the square of the coupling and decrease with the squares of the interaction duration, the energy 
splitting between the parallel energies and the tilt (chirp) rate. Physical examples of this model can 
be found in coherent atomic excitation and optical shielding in cold atomic collisions. 

PACS numbers: 32.80.-t, 32.80.Bx, 33.80.-b, 34.50.-s, 33.80.Bc, 32.80.Qk 



I. INTRODUCTION 

The famous Landau-Zener (LZ) model [l| is the most 
popular tool for estimating the transition probability be- 
tween two states whose energies cross in time. This model 
assumes a constant interaction of infinite duration and 
linear energies. Owing to some mathematical subtleties, 
the LZ model often provides more accurate results than 
anticipated (given the rather simple time dependence of 
the LZ Hamiltonian) when applied to real physical sys- 
tems with more sophisticated time dependences. The 
popularity of the LZ model is further motivated by the 
extreme simplicity of the transition probability. 

The LZ model has been extended to three and more 
levels by a number of authors. There are two main types 
of generalizations: single-crossing (bow-tie) models and 
multiple-crossings grid models. 

In the bow-tie models all state energies cross at the 
same instant of time. Carroll and Hioe have solved the 
three-state bow-tie model in a special symmetric case 0] 
and in the general case An extension of this model 
to N states has been suggested 0, Q and then rigorously 
derived by Ostrovsky and Nakamura 0. A further ex- 
tension, wherein one of the levels is split into two parallel 
levels, has been suggested Q and derived [1] by Demkov 
and Ostrovsky. An example of a bow-tie system occurs 
when a sequentially coupled quantum ladder of states in 
an atom or a molecule is driven by chirped laser pulses 
0; an adiabatic sweep of frequency through resonance 
will transfer all population from the lowest to the highest 
energy state, for either sign of the chirp. A bow-tie type 
of linkage can also occur in a rf-pulse controlled Bose- 
Einstein condensate output coupler [l(| [H| • Yet another 
example is the coupling pattern of Rydberg sublevels in 
a magnetic field 

In the multiple-crossings grid models the energy dia- 
gram consists of a grid of crossings formed by two man- 
ifolds of rectilinear parallel diabatic energies that cross 
each other. In the Demkov-Osherov (DO) model [H|,[l3j], 
a single tilted energy crosses a set of N parallel energies. 
This model has been generalized to the case when the 



single tilted energy is replaced by a set of M parallel en- 
ergies, which cross the other set of N parallel energies 
[14l - [l6j . The special case when M and N are infinite (so 
that the grid of crossings is periodic) has also been solved 
17]. Effects of level degeneracies and quasi-degeneracies 
have been studied by Yurovsky and Ben-Reuven [3, [lj| . 

In the most general case of an asymmetric linear 
Hamiltonian, H(i) = A+Bi, where B is diagonal, the gen- 
eral solution has not been derived yet, but exact results 
for some "survival" probabilities have been conjectured 
d and derived [HQ. 

A variety of physical systems provide examples of mul- 
tiple level crossings. Amongst them we mention ladder 
climbing of atomic and molecular states by chirped laser 
pulses [HI , optical shielding in cold atomic collisions (25r - 
29j . optical centrifuge for molecul es 13011 . Stark-chirped 
rapid adiabatic passage (SCRAP) |3lT[33| . and creation 
of entaglement in many-particle systems [34j , In fact, 
physical situations where multiple level crossings play a 
role have been discussed already in 1960's, when the har- 
poon model for reactive scattering was considered (35[ . If 
one adds the vibrational states to the scattering picture, 
then one obtains a multilevel crossing model, which has 
been studied using the LZ model [3o| . 

A general feature of all soluble multilevel models is that 
the transition probabilities P m ^ n between states ip m and 
if) n are given by very simple expressions, as in the origi- 
nal LZ model, although the derivations are usually quite 
involved. In the grid models in particular, such as the 
DO model, the exact probabilities P m -¥n have the same 
form — products of LZ probabilities for transition or no- 
transition applied at the relevant crossings — as what 
would be obtained by naive multiplication of LZ proba- 
bilities while moving across the grid of crossings from ip m 
to tp n , without accounting for phases and interferences. 

A very interesting feature of all grid models is that 
counterintuitive transitions, for which the level crossings 
appear in a "wrong" order, are forbidden. In the three- 
state example in Fig. [T]such is the transition ip^ — > tpi. 
In the adiabatic limit, the inhibition of such transitions 
is easily understood: for the ip3 — ¥ ipi transition an adia- 
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FIG. 1: (Color online) Diabatic and adiabatic energies in 
the three-state crossing model studied in this paper. We 
are concerned primarily with the (counterintuitive) transition 

^3 ~> tpl- 



analytical estimates for it. Our approach involves trans- 
formation to the adiabatic basis where the evolution is 
represented as a sequence of instantaneous two-state LZ 
transitions at each crossing and adiabatic evolution else- 
where. This approach allows us to derive the transition 
probabilities between each pair of diabatic states, includ- 
ing the probability for the counterintuitive transition. 

The problem of counterintuitive transitions, besides 
quite interesting by itself, has interesting physical im- 
plications, which are discussed in some detail in Sec. I VII 
Amongst them, we mention the problem of saturated op- 
tical schiclding with near-resonant light, which plays an 
important role in cold atomic collisions [25 29]. 

This paper is organized as follows. We define the prob- 
lem in Sec. Hi] and the propagator is derived in Sec. IIIII 
in the general case. The transition probabilities in the fi- 
nite DO model are derived and compared with numerical 
results in Sec. IIVI The time-dependent probabilities in 
the original (infinite) DO model are presented in Sec. [V] 
Some physical examples of counterintuitive transitions 
are discussed in Sec. IVII The conclusions are summa- 
rized in Sec. IVII1 



batic path requires the level crossing between ipz and 
to occur before the crossing between V2 and this is 
not the case here and hence there is no adiabatic path 
linking 1^3 to tpi . In the general case of nonadiabatic evo- 
lution, however, the inhibition of the tp3 ~~ ^ "01 transition 
is not so obvious because the concerned final state ipi 
acquires some nonzero transient population during the 
interaction. Yet in the end this population vanishes, and 
this is an exact result. A similar conclusion applies to 
more general level-crossing models [2l| [23| as well. 

We have verified with numerical simulations that all 
ingredients of the multilevel LZ models are essential for 
this feature: linear energies, constant interactions and 
infinite duration. Nonlinear energies, pulsed interactions 
or finite interaction duration can each lead to nonzero 
probability for counterintuitive transitions. 

Yurovsky et al. [361 ] have studied analytically and nu- 
merically the counterintuitive transition probability in 
two variations of the DO model: (i) with a finite interac- 
tion duration and (ii) with a piecewise-linear sloped po- 
tential. They have used a perturbative approach assum- 
ing a quasidegenerate band of parallel energies and have 
found nonzero probabilities for counterintuitive transi- 
tions in both models. 

In this paper, we derive analytically the probability 
for a counterintuitive transition in the simplest case of 
the DO model involving three states, with two parallel 
and one slanted energy, as shown in Fig. [TJ We assume 
a finite interaction duration, but our approach does not 
use the quasidegeneracy assumption of Yurovsky et al. 

therefore our results are more general, as far as the 
three-state case is concerned. The purpose of this work 
is not just to show that the probability for a counterintu- 
itive transition is nonzero but rather to derive accurate 



II. DEFINITION OF THE PROBLEM 

A. The system 

The probability amplitudes of the three-state system 
C(t) = [Ci(t), C<z(t), C3(t)] T with the energies shown in 
Fig. [1] satisfy the Schrodinger equation (K = 1), 

iC(t) = H(t)C(i), (1) 

where the overdot denotes d/dt. The Hamiltonian in the 
usual rotating-wave approximation is given by [37| 
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As the Hamiltonian ^ shows, there is no direct coupling 
between states ipi and ^3 but each of them is coupled to 
state ip2 ■ Without loss of generality the constant cou- 
plings f2i2 and Q23 will be assumed real and positive. 
Both couplings are supposed to have the same finite du- 
ration, being turned on at time i$ and turned off at time 
tf. In the original DO model the couplings last from —00 
to +00. Furthermore, the energy splitting parameter 5 
and the slope (5 of the energy of state ip2 are assumed 
positive too, 

6 > 0, > 0. (3) 

Given these assumptions, the crossing between the di- 
abatic energies of states ipi and ip2, occuring at time 
t- = —t precedes the crossing between states ip2 and 
tp3, occuring at time t+ = r, where 
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Therefore, the transition ip\ — ¥ -03 is intuitive, while the 
opposite transition ip 3 — > ipi is counterintuitive. In the 
adiabatic limit, the transition probability from tp\ to 1^3 
is Pi_>3 = I, whereas that from ip 3 to ipi is P^i = 0. 



B. The Demkov-Osherov model 

In the DO model the transition probabilities P m ->« 
from state ip m at t — > —00 to state ip n at t — > +00 are 
given exactly by products of two-state single-crossing LZ 
probabilities, as follows 

Pl-h =p~, P\^2 = q-p+, Pi^3 = q-q+, 

Pz-yi = q-, P 2 ^2=P-P+, P 2 ^3=P-q+, (5) 
P3-H = 0, P 3 ^2 = q+, Ps^3 = P+, 



where 



where 



P± 



-2ira^ 



q± 



1 



(6) 



i.e. g± is the transition probability and p± is the proba- 
bility of no transition at the crossing t± , with 



(7) 



These simple results coincide with what would be ex- 
pected naively, by treating the crossings independently, 
no matter how close they are to each other, and multi- 
plying LZ probabilities. In particular, if the system is 
initially in state ^3, the transition probability to state 
ipi is exactly zero at t — > +00, = 0, which means 

that the counterintuitive transition ip 3 — > ipi is forbid- 
den. This zero probability is rather unexpected because 
state ipi acquires some nonzero population during the in- 
teraction. However, it vanishes at t — > +00 for any set 
of parameters, irrespective of whether the interaction is 
adiabatic or not. This property is unique for the DO 
model and it depends crucially on any of its features: 
infinite coupling durations, constant couplings, constant 
energies of states ip% and V>3 and linear energy of state ijj 2 - 
The goal of the present paper is to estimate the proba- 
bility for counterintuive transitions in the case of finite 
coupling duration. 



III. TRANSITION MATRIX 

A. Eigenvalues and eigenstates 

We need the eigenvalues and the eigenstates (the adi- 
abatic states) of H(t). The eigenvalues read [13] 



1 2 1 

Ai = 3 SCOS 3^' 

12 1 
A 2 = --a- -scos-(f + 7r), 

12 1 
A 3 = --a - -scos-(6> - tt), 



(8a) 
(8b) 
(8c) 



a = -fit, 

b = -(s 2 + ci 2 12 + n 2 23 ), 
c = s(n 2 12 -n 2 23 + spt), 



cos 9 



s = \/ a 2 - 36, 

2a 3 - 9a6 + 27c 



2s 3 



(9a) 
(9b) 
(9c) 

(9d) 
(9e) 



The eigenstates are given by ipk = [fik, hk, hk] T , with 



hk = tUa 2 ^ 2 ), 
J^k 

hk = ^r^2z{\k + 5), 



(10a) 
(10b) 
(10c) 



where Nk are normalization factors (k — 1,2,3). The 
asymptotic behaviors at large times of the eigenvalues 
and the eigenstates are presented in Appendix [X] 



B. Adiabatic basis 

The transformation linking the diabatic amplitudes 
C(t) and the adiabatic amplitudes A(t) is given by 



C(t) = F(t)A(t), 



(11) 



where A(t) = [A^t), A 2 (t), A 3 {t)] T and F(t) is an orthog- 
onal rotation matrix [F _1 (i) = F T (i)] whose columns are 
the eigenvectors (fTT 



F(t) 



fn(t) Aa(t) f 13 (t) 
hi(t) f 22 (t) f 23 (t) 

hi(t) f 32 (t) / 33 (t) 



(12) 



The Schrbdinger equation in the adiabatic basis reads 

iA(t) = H A (t)A(t), (13) 
with H A (t) = F T (t)H(t)F(t) - iF T (t)F(i), or 



Ha(*) = 



Ai -ivu -ivi 3 
-iv 2X \ 2 -iv 2Z 
-ivz\ ~iv 32 A 3 



(14) 



where the nonadiabatic coupling between the adiabatic 
states <pk(t) and (pi(t) is 



v k i{i) = (ifik(t) I tpi(t)) = -Hk{t) 



(15) 



We use the fact that the transition times in the adi- 
abatic basis are shorter than in the diabatic basis [381 ) . 
This is so because while the asymptotic behaviors of the 
adiabatic energies at large times (|A2[) are approximately 
the same as the asymptotics of the diabatic energies, the 
couplings Vki in the adiabatic basis (fT5")) vanish as t~ 2 [see 



Eqs. (|A2[) ], in contrast to the constant couplings fli 2 and 
^23 in the diabatic basis. The difference in the transition 
times is illustrated in Fig. [51 where the oscillations in the 
populations of the adiabatic states vanish much faster. 
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FIG. 2: (Color online) Time evolutions of the populations of 
the diabatic and adiabatic states for the original DO model 
(ti — — oo). The system starts in state tpi. The interaction 
parameters are fii2 = f^23 = j3? , 5 — ft? . 



C. Evolution matrix in the adiabatic basis 

Our method is based on two simplifying assumptions. 
First, we assume that appreciable transitions take place 
only between neighboring adiabatic states, (fi(t) ■<-> <f2(t) 
and (p2(t) -n- if3(t), but not between states p\(t) and 
(y9 3 (t), because the energies of the latter pair are split by 
the largest gap. Second, we assume that the nonadiabatic 
transitions occur instantly at the corresponding avoided 
crossings and the evolution is adiabatic elsewhere. This 
allows us to obtain the propagator in the adiabatic basis 
by multiplying five simple transition matrices describing 
LZ transitions or adiabatic evolution. 

The adiabatic evolution matrix U A (tf,ti) is most con- 
veniently determined in the adiabatic interaction repre- 
sentation, where the diagonal elements of H^i) are nul- 
lified. The transformation reads 



A(t) = M(t)B(i), 



KMo) = 



-»Ai(t,t ) 








; -tA 3 (t,to) 







,-iA 3 (t,t ) 



A fc (Mo) = / Xk(t')dt\ 
A k i(t,t ) = A fc (t,t )- M(t,t ), 



(17) 



(18a) 
(18b) 



and to is an arbitrary fixed time. The Schrodinger equa- 
tion in this basis reads 



with 



H B (i) = -i 



u 21 e 



*B(t) = H fl (t)B(t) 


iA 2 i(Mo) 



(19) 



j/ 12e *Ai2(*,to) jy 13e iAi 3 (t,to) 







^23 e 



iA 23 (t,*o) 



I/31 giA 3 l(t,to) j, g«A 3 2(t,to) 







In this basis, the evolution matrix for adiabatic evolution 
is given by the identity matrix. 

The LZ transitions at the crossings at ±r are described 
by the transition matrices 



Ulz(-t) 
Ulz(t) 



1 

y/qle-**- 
JpZ 





fqZe 



P- 



q+e -VPT 

y/q+e 1 ^ 
1 



(21a) 
(21b) 



where p± and q± are given by Eqs. ([5]) and 



b ± = argT(l - ia±) + - + a 2 ± (lna| - l) 



(22) 



(16) 



with a± given by Eqs. (0. The LZ phases 4>± do not 
depend on time, unlike the dynamical phases (|18al) . 
The propagator in the adiabatic basis reads 

U A {tf,ti) = M(tf,T)U LZ (T)M(T,-T)U LZ (-T)M(-T,U), 

or 



U A (t f ,U) = 



fq^ e - i( t>+-i^x{tf,U) _ /p + q_ e -#- -» A i (tf i T ) -» A 2 (r,U ) ^/ p_p + e - iA i(tf > T )-' iA 2 (T,-T)-iA 3 (-T,ti) 
pp—- e -iA 1 (T,t i )-iA 2 {t s ,T) fq— g— e i(0 + -^-)-iA 3 (t/,ti) _^ / p_q + e ief> + -iA 2 {tf,-T)-iA3(-T,t i ) 







/p_e 



-iA 2 (-T,ii)— {As (t/,-r) 



ic/> - -i A 3 (tf,ti) 



.(23) 

In the special case when = — T, tf = T and f2i2 = ^23 = many expressions simplify, as shown in Appendix |B| 
Then a + = a_ = a, p + = p- = p, q + = q- = q = I — p, and Eq. f|23[) reduces to 



U A (T,-T) 



/q e -i<l>-iA x {T-T) ~^/pq e -i0-iA 12 (T,T) 
^ e iA 3 (T,-r)-<A 2 (T,r) g e - 

^/p e iA 2 (T,r)-iA 3 (T,-T) e #+iAi (T,-T) 



o i0+iAi 2 (T,r) 



(24) 
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D. Evolution matrix in the diabatic basis 

The propagator in the diabatic basis can be obtained 
by using the transformation it reads 



0.4 



U(t f ,t i ) = F(t f )U A (t f ,t i )F T (t i ). 



(25) 



We shall use this relation to derive the transition proba- 
bilities in the finite and original DO models below. 



IV. TRANSITION PROBABILITIES IN THE 
FINITE DEMKOV-OSHEROV MODEL 

A. The propagator 

In order to obtain simpler formulas for the proba- 
bilities we assume that ti = — T, tf = T, although 
our approach is not limited to these restrictions. The 
transition probability from state to ip n is given by 



P m _j. n = \U nm (T,-T)\ , where 



3 



U nm (T, —T) = \ f nk (T)U£(T, -T)f ml (-T). (26) 



k,l = l 



Using this relation one can calculate the transition prob- 
ability between any two states of the system. We pay 
special attention to the probability for counterintuitive 
transitions Pa-^i, which is zero in the original DO model. 



B. Counterintuitive transition 

In the special case of equal couplings, f2i2 = ^23 = ^, 
we find from Eqs. ((Ml), (Ell) and (|B3)l that the transition 
probability F 3 ^i = \U 13 (T,-T)\ 2 reads 

Pb-H = \pfn+qfu-2^f 11 f 12 cos[A 12 (T,T) + ( f>} 
+2^/12/13 cos[A 2 (T, t) - A 3 (T, -t)] 
+2^f 11 fi3CO S [A 1 (T,-T) + cP}\ 2 . (27) 



It can be written as 



P 



3^1 



P 



3-s-l 



P. 



3-s-l, 



(28) 



where P 3 ^,\ is the average probability and P 3 ^i is the os- 
cillating part. By using the asymptotic expansions (|A3p 
for f mn in Appendix [X] and keeping the leading terms in 
the expansion over 1 jT we find 



P 



3-^1 



(pfn + qfh) 2 + %/n + pfh)fu + zpqf&9& 

,(29b) 



n 2 (n 2 P + A5 2 q ) n 2 (n A P - 86 4 q ) 

"T" f» r-Q /-iQm 1 ? "T" 



2S 2 /3 2 T 2 



2<5 3 /3 3 T 3 



numerical 



average analytical 

analytical 




2 4 6 8 10 

Interaction Duration T (units of p~' /2 ) 

FIG. 3: (Color online) Probability for counterintuitive tran- 
sition P3_).i plotted against the coupling duration T for inter- 
action parameters Q12 = Q23 = , S = (3?. The full line 
shows numerical results, the thin dashed curve the analytical 
approximation (|270 and the thick dashed curve the average 
probability ((29}. 



CI 2 



p 2 T 2 



{4qcos 2 [Ai(T,-T) + , 



77$1 2 

_^_ cos 2 [A 2 (T,r)~ A 3 (T,-r)] 
AJpqtt 



Pv 



cob[Ai(T,-T) + 0] 

x cos[A 2 (T,r)-A 3 (r,-r)]} + ... (30) 

For unequal couplings (f2i2 ^ ^23), the expansion over 
1/T of the average probability reads 

n 2 2 nh(p- + p+) + 45 2 (n 2 12 q + + n 2 3 q-) 



p 



3->l 



45 2 /3 2 T 2 

n 2 2 n 2 3 (n 2 2P _ + n 2 23P+ ) - 8S%n 2 2 q+ + n 2 23 q_) 



4<5 3 /3 3 T 3 



(31) 



The part P3_>i is too cumbersome to be presented here. 

In the near-adiabatic regime (a ^> 1, p <C 1, q fa 1, 
and <p < 1), we find from Eq. $ZT§ that 



3-S-l 



Ail 2 
^T 2 



cos 2 [A 1 (T,-T) + . 



(32) 



In the weak- coupling limit (a <C 1, p ~ 1, 9<C 1, and 
w 7t/4), we find from Eq. ([57]) that 

P ^ ~ ^a^^M^-Aa^-T)]. ( 33 ) 

Figure [3] shows the counterintuitive transition proba- 
bility P3_n against the coupling duration T. The proba- 
bility decreases in an oscillatory manner, as predicted by 
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our results. The analytical approximation ([27)1 describes 
very accurately both the phase and the amplitude of the 
oscillations. The approximation (|29|) describes very ac- 



curately also the average probability P3-H. 

Fig. [4] displays the counterintuitive transition prob- 
ability -P3-).i as a function of the energy separation pa- 
rameter 5. As with T, the probability decreases in an 
oscillatory manner. The analytical approximations are 
seen again to describe the probability very accurately. 



C. Other transition probabilities 

By using Eq. (|26[) we can find all transition proba- 
bilities P m _,„ = \U nm (T,-T)\ 2 (m,n = 1,2,3). For the 
sake of simplicity we assume again that U = —T, tf = T 
and O12 = f^23 = ^, although our approach applies to 
the general non-symmetric case as well. Using Eq. (f2"4")l 
we find the average probabilities expanded to the lowest 
order of 1/T, 



Pi- 



tt 2 



Pi- 



Px- 



P2- 



P 



2^2 



P2- 



P 



3^1 



P + 
pq 

I 2 - 

q + 
p 2 - 

pq + 

2Q 2 

(3 2 T 2 



f3 2 T 2 
tt 2 



[^(^-2p) + l-2p-^],(34a) 



/3 2 T 2 

2tt 2 

P 2 T 2 
tt 2 



{x 2 q 2 + l-6p + 7p 2 ), 
[x 2 (p-q 2 )+3pq-q], 



p2 T 2 
2fl 2 

~ p2 T 2 

n 2 



(p 2 +Ap-3-x 2 q 2 ), 
(1+p-Ap 2 ), 
(x 2 q 2 + l-6p + 7p 2 ) 



P2 T 2 

(x 2 p + q), 



P3- 



tt 2 



P 



3^3 



q 

p- 



P2 T 2 

n 2 

p 2 T 2 



(p 2 +Ap-3-x 2 q 2 ), 



(34b) 
(34c) 
(34d) 
(34e) 
(34f) 
(34g) 
(34h) 



x 2 (q 2 -2p) + l-2p-p%(Mi) 



where x 2 = £l 2 /4c> 2 . All these probabilities have the 
correct DO limits ^ for T — > 00. Note that Pi_>i ~ 

Ps^3, Pl^2 ~ P2^3, and P 2 ^i ~ P3-S-2- 



V. TIME EVOLUTION IN THE ORIGINAL 
DEMKOV-OSHEROV MODEL 

In the original DO model the time-dependent tran- 
sition probability from state ip m to i\) n is given by 
P m _> n (t) = I [/„,„(£, -00) | 2 . We find from Eq. for 
U = —00, tf =t that 

3 

U nm (t, -00) = ^ fnk(t)Uki(t, -oo)/ mt (-oo), (35) 
fe,/=i 



numerical 



average analytical 

- analytical 




2 4 6 

Splitting 8 (units of P' /2 ) 



FIG. 4: (Color online) Probability for counterintuitive tran- 
sition P3_>i plotted against the energy splitting parameter 
8. The other interaction parameters are Q12 = Q23 = /3 5 , 
T — 5/3~2. The full line shows numerical results, the thin 
dashed curve the analytical approximation (|27p and the thick 
dashed curve the average probability (|29|) . 



and we take into account that 



F(-oo) = 



0-10 

1 

1 



(36) 



Counterintuitive transition 



By using Eq. (|35[) we find that the probability of coun- 
terintuitive transition P3->i after the crossings reads 



P 3 ^i(i) = lv^7ii(t)+VSf/i2(*)e 



i[Al2 (t,r)+<p + 



p 3 ^l(t) + p 3 ^l(t) 



where 



P 3 ^i(t) = < ?+ / 2 1 +p+/ 2 2 , 

^12 {^ 2 q+ + tt 2 ^ 



A5 2 p 2 t 2 



(37b) 

(38a) 
(38b) 



P 3 ^i(<) = 2vpT^/ii/i2 coB[Aia(t,T) + 0+] (39a) 

^12^23 



5f3 2 t 2 



■^p+q+ cos[Ai 2 (i, r) + <f>+\. (39b) 



Figure [5] shows the time evolution of the probability 
for counterintuitive transition Ps_>i for three values of 
the couplings O12 = = O. For these values of 57 
there are almost no oscillations visible, because p+ <C 1 
in Eq. (|39bl) . The probability p3_ > i(t) decreases with 
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time, as predicted by Eq. (|38|) . and vanishes for large 
times towards the DO result P 3 ^i(oo) = 0. As predicted, 
the transition probability P 3 _ > i(i) increases with the cou- 
plings. The analytical approximation (|38l) describes very 



0.20 



accurately the average probability P^i(t). 



B. Other transition probabilities 

By using Eq. and the large-i asymptotic expan- 
sions in Appendix [XJ we find the leading terms of all 
average transition probabilities P m ^„(t) after the cross- 
ings for equal couplings (fii2 = ^23 = 

Q 2 





+l(*) 


Pl- 




Pl- 


+s(t) 


P2- 


+i(i) 


P 2 - 




P 2 - 


+s(t) 


P3- 




P 3 ^2(t) 


P3- 


,3(0 



pn 2 

n 2 
f p 2 t 2 1 

n 2 



1 1 (3 2 t 2 LH * * 

/ 2 

1+05*2 (P "9- 

02 + t£( 1 -3p 2 ) 

fl 2 r , 



,2„2\ 



,2„21 



/3 2 i 2 



(xr 2 p + g) , 



n 2 

/3 2 t 2 ' 

r; 2 

^ 2 t 2 



(40b) 
,(40c) 
(40d) 
(40e) 
(40f) 
(40g) 
(40h) 
(40i) 



where, as before, x 2 — Vt 2 /A8 2 . All these probabilities 
have the correct DO limits (J5J for t — > 00. It is interesting 
to note that the counterintuitive-transition probability in 
the original DO model Q40g[ ) is one-half of the one in the 
finite DO model {Mg}. 



VI. EXPERIMENTAL IMPLEMENTATIONS 

This three-state DO model discussed here can be re- 
alized in several physical systems. One example is when 
a chirped laser pulse couples an initially populated state 
to a manifold of two levels simultaneously [9| , the energy 
diagram for which is identical to the one in Fig. [TJ with 
state ip2 being the initial state. Adiabatic evolution can 
produce in this system a very selective excitation even if 
the Fourier bandwidth of the laser pulse is larger than 
the level spacing within the manifold. Indeed, for red- 
to-blue chirp (/? > 0) the system will follow from left 
to right the lowest adiabatic energy, which links the ini- 
tial state ip2 to the lowest diabatic (unperturbed) state 
-01 of the manifold. In contrast, for blue-to-red chirp 




analytical 
numerical 



10 20 30 

Time (units of p~' /z ) 



40 



FIG. 5: (Color online) Time evolution of the probability for 
counterintuitive transition P3_>i in the original DO model 
(ti — — 00) for 5 = /9* and three values of the couplings, 

Q12 = ^23 = /? 5 , 3/?2, and 10/3 3 (denoted near the respec- 
tive curves) . The solid curves show numerical results and the 
dots show the analytical approximation (|38[) for the average 
probability. 



(/? < 0) the system would follow (from right to left in 
Fig. [T]) the highest adiabatic energy, linking state tp2 
to the highest diabatic state tpz- Therefore, the chirp 
sign alone determines if the population is directed to- 
wards the lowest or the highest state of the manifold. 
This excitation scheme has been demonstrated experi- 
mentally by Warren and co-workers [39j on the 3s-3p 
transition in sodium. Red-to-blue chirped picosecond 
pulses populated predominantly the lower fine-structure 

1/2; 



level 3p 2 Pi/2, while blue-to-red chirped pulses placed the 
population onto the upper fine-structure level 2>p 2 P^,/2- 
Counterintuitive transitions can be demonstrated in this 
system by using the same setup as described above, but 
starting in one of the fine-structure levels. 

Another example of the model discussed here is Stark- 
chirped rapid adiabatic passage (SCRAP) [3ll - [33j where 
level crossings are created by inducing ac Stark shifts of 
the energy levels by a strong off-resonant laser pulse. A 
level diagram similar to the one in Fig. [T]is found in the 
three-state version of SCRAP [33|. 

It is also worth noting a similarity between the fi- 
nite DO model presented here and multiple-ground-state 
models of optical shiel ding of cold collisions in magneto- 
optical atom traps |25H29j . Optical shielding techniques 
were originally developed to prevent laser cooled atoms 
from escaping the trap by making the collisions between 
the atoms elastic [26|. The experiments show that the 
shielding process saturates, in contrast to the expected 
complete shielding, when the intensity of the shielding 
field is increased [2a]. The reason for the saturation 
is still an open problem and may have a connection 
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to the counterintuitive transitions in the corresponding 
multiple-ground-state level scheme of the quasimolecule, 
which in its simple form resembles the present model. 



2001-00294. JP acknowledges support from the Mag- 
nus Ehrnrooth Foundation and thanks for the hospitality 
during his visit to Sofia University. 



VII. CONCLUSIONS 

Wc have calculated analytically the probability for the 
counterintuitive transition — > ipi in a three-state sys- 
tem with crossing energies. We have performed the 
derivation in the adiabatic basis by assuming instanta- 
neous transitions at the level crossings and adiabatic evo- 
lution elsewhere. The counterintuitive transition proba- 
bility P 3 _>i is nonzero at finite times, whereas it vanishes 
for infinite coupling duration, thus recovering the result 
in the DO model. A very good agreement with the nu- 
merical results is found. This approach has been used to 
derive the other transition probabilities within the three- 
state system too. 

Our results suggest that at large T the probability for 
a counterintuitive transition is proportional to the factor 
O 2 / (S 2 f3 2 T 2 ) . The decrease with the interaction duration 
T is expected as for T — > oo the probability should vanish 
(DO model). The increase with the coupling O is also an- 
ticipated since larger interaction is expected to increase 
weak transitions. The decrease vs the energy splitting 5 
between the concerned states ipi and ip3 is natural too. 
The decrease vs the slope f3 is more subtle. Since (3T 
is the diabatic energy of state ip%, and since states tpi 
and ips interact with each other via ip2 , the effective cou- 
pling duration for the transition ip 3 — > tpi depends on 
the slope /3: the larger the slope the smaller the duration 
and therefore the smaller the transition probability. The 
other ("intuitive") transition probabilities exhibit simi- 
lar dependences of the finite-duration corrections on the 
interaction parameters. 

We have described physical examples of coherent 
atomic excitation where the described counterintuitive 
transition can be observed. The model studied here re- 
sembles also the simple form of the multiple-ground-state 
quasimolecule model used to study optical shielding in 
magneto-optical atom traps [25L [26j j . The saturation of 
optical shielding in high laser intensities has been earlier 
attributed to off-resonant excitation to attractive states, 
counterintuitive transitions or other processes including 
multiple partial waves [25U27I ] . The results here show 
that the probability for counterintuitive transitions is 
clearly non-negligible and indicate that they can also play 
a role in the saturation of optical shielding. 
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Appendix A: Asymptotics of the eigenvalues and the 
eigenstates 



Here we present the asymptotics of the eigenvalues 
([5]) at large positive time t > 0, i.e. for j3t ^> S and 
pt > O. Then a = -fit, b = -(5 2 
c = 8 2 pt [l + (0 2 2 - 2 3 ) /(Spt)] , and 



O 2 
"12 



^23)' 



cos 9 



Pt 
1 - 



3(s 2 + o 2 12 + o 2 3 ) 



2/3t 



27S 2 



2/3 2 t 2 ' 
Pt ' 



0\ 2 



^23 



The eigenvalues have the asymptotics 
Ai - pt 

a 2 ~ s- 



pt 



^23 

Pt ' 



A, 



O 2 

Pt 



Hence we find from Eqs. (fT0| that 



<Pi 



9'2 



¥>3 



1 - 



flu 

pt 

"12^23 



0,12 + ^23 ^23 

2p 2 t 2 ' Pt 



2Spt ' 



^23 2 23 (0 2 2 +45 2 ) 

pt ' 85 2 P 2 t 2 



(Ala) 
(Alb) 
(Ale) 

(A2a) 
(A2b) 
(A2c) 

(A3a) 

n T 

(43b) 



AS 2 



8S 2 P 2 t 2 



O12 

pt ' 



^12^23 

28 pt 



(A3c) 



Appendix B: Symmetries of the eigenvalues and the 
eigenstates 

The eigenvalues and the eigenstates simplify when = 
— T, tf = T and O12 — ^23 = O. Then a + = a_ = a, 

p+ = p- = p, q+ = q- = q = l - p, a{T) = -pT, 

b(T) = -5 2 - 20 2 , c(T) = S 2 PT, and 



s(T) = ^P 2 T 2 + 3S 2 + 6O 2 



cos 0{T) 



PT 



{P 2 T 2 - 9S 2 + 90 2 ). 



(Bla) 
(Bib) 



Hence 6{-T) = n - 9(T) and therefore A x (-T) = 
-A 3 (T), A 2 (-T) = -A 2 (T), A 3 (-T) = -A X (T), and 

A 2 (-r,-T) - A 2 (r,T) = -A 2 (T,r), (B2a) 

A 2 (T,-T) = A 2 (r,-r)=0, (B2b) 

Ai(-r,-T) = A 3 (r,T) = -A 3 (T,r) ) (B2c) 

A 3 (-r, —T) = A, (r, T) = -Ai (T, r) . (B2d) 
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The transformation matrix at — T is given by 



F(-T) = 



-f 33 (T) -/ 32 (T) -/ 3 i(T) 
/ 23 (T) / a2 (T) / 21 (T) 
-/ia(T) -/ 12 (T) -f n (T) 



With these relations taken into account the propagator 
([25]l reduces to Eq. (jSJ). 



(B3) 
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